SUPERNOVA SN 1987A REVISITED AS A MAJOR PRODUCTION 

SITE FOR r-PROCESS ELEMENTS 
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Abstract 

The origin of nucleosynthesis products of rapid neutron capture reactions (the r-process) is a long- 
standing astrophysical problem. Recent analyses of elemental abundances for extremely metal-poor stars 
shed light on the elemental abundances of individual supernovac. Comparison of the abundance distri- 
butions of some extremely metal-poor stars with those of the best-observed supernova SN 1987A clearly 
indicates that the overabundances of barium and strontium found in SN 1987A that have been ascribed 
to the slow neutron capture process must be results of r-process nucleosynthesis. The mass of freshly 
synthesized barium in SN 1987A is estimated to be 6 x 1O~ 6 M0 based on the observed surface abundance 
and detailed hydrodynamical models for this supernova. These new findings lead to the conclusion that 
20 Mq stars, one of which is the progenitor star of SN 1987A, are the predominant production sites for 
r-process elements in the Galaxy and the r-process element donors for notable ncutron-capture-rich giant 
stars, CS22892-052 and CS31082-001. 
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1. INTRODUCTION 

The existence of astrophysical r-process elements (ele- 
ments synthesized by rapid neutron capture process) in 
old metal-poor stars has indicated that these elements 
are synthesized in the preceding supernovae (SNe) orig- 
inating in massive stars (Truran 1981). The elemental 
abundance determination of numerous Galactic halo stars 
with extremely low metal abundances has opened the door 
to much more detailed investigations of r-process produc- 
tion from SNe in the early Galaxy (McWilliam ct al. 1995; 
Ryan, Norris, & Beers 1996). However, which SNe pro- 
duce the r-process elements has not yet been agreed upon 
in spite of many authors' attempts (Mathews, Bazan, 
& Cowan 1992; Ishimaru & Wanajo 1999; Tsujimoto, 
Shigeyama, & Yoshii 2000; Qian 2001; Fields, Truran, & 
Cowan 2001). 

Direct information about the production site for re- 
process elements should come from abundance determi- 
nation of the elements in the ejecta of individual SNe. 
Absorption lines due to barium (Ba) and strontium (Sr), 
~10% of which are of r-process origin in the solar ele- 
mental abundance, were detected in the spectra of SN 
1987A (Williams 1987). A remarkable feature found by 
Mazzali, Lucy, & Butler (1992) was that the shape of 
the Ba absorption line indicated lack of this element at 
the very surface of the ejecta. Thus the observed Ba 
must have been synthesized inside the star and did not 



exist in the interstellar matter (ISM) from which the 
star formed. The abundances of Ba and Sr were de- 
rived by taking into account non-LTE effects. The result 
showed that these elements are overabundant in compar- 
ison with the abundance of the Large Magellanic Cloud 
(Hofhch 1988; Mazzali, Lucy, & Butler 1992; Mazzali & 
Chugai 1995). The overabundances of Ba and Sr in SN 
1987A have led a few authors to argue that these ele- 
ments were synthesized by a slow neutron capture process 
(weak s-process) operating in the core helium (He) burn- 
ing stage of the progenitor star (e.g., Prantzos, Arnould, 
& Casse 1988). However, the derived abundance ratio 
of these two elements in SN 1987A was larger than the 
corresponding solar ratio, i.e., Ba/Sr~2.5(Ba/Sr) Q . This 
value is clearly inconsistent with the results of weak s- 
process calculations (Prantzos ct al. 1988) which pre- 
dict a ratio significantly smaller than the solar ratio such 
as O.l(Ba/Sr) <Ba/Sr<O.6(Ba/Sr) (Fig. 1). Another 
fact casts doubt on s-process origin: other SNe II-P 
such as SN 1985P, SN1990E, and SN 1990H did not 
show any sign of the overabundance of Ba (Chalabaev & 
Cristiani 1987; Mazzali & Chugai 1995), though s-process 
calculations predict the overabundance of Ba irrespective 
of the He core mass (Prantzos et al. 1988). Together, these 
facts strongly suggest that these elements in SN 1987A are 
not of s-process origin. 

On the other hand, the derived Ba/Sr ratio can be 
reconciled with the ratios of the same elements in two 
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Fig. 1. Comparison of the elemental abundance distribu- 
tion in the ejecta of SN 1987A with that of the star 
CS22892-052. For Sr and Ba, the abundance distribu- 
tion of another neutron-rich giant star CS31082-001 is also 
plotted together with the solid and dashed lines showing 
the abundance distributions theoretically anticipated from 
the r-process origin and s-process nucleosynthesis (Prantzos, 
Arnould, & Casse 1988), respectively. These lines are 
drawn from the [Ba/Fe] corresponding to the ratio of 
CS22892-052. For lighter elements between Mg and Ni, 
the abundance curves predicted by a recent nucleosynthesis 
calculation for 20AfQ (solid line) and 13Mq (dashed line) 
stellar models (Umeda & Nomoto 2001) are also shown. 



extremely metal-poor stars CS22892-052 and CS31082- 
001 (Fig. 1). Since these elements in such metal- 
poor stars must be results of r-process nucleosynthesis 
(McWilliam 1998), the derived ratio suggests that these 
elements in SN 1987A might also be of r-process origin. If 
this is the case, it is expected that the observed Ba (Sr) is 
synthesized during explosion in the deepest layers of the 
ejecta where the matter is exposed to an intense flux of 
neutrons. At the same place, a radioactive element nickel 
56 ( 56 Ni), which is one of the major seed elements for 
r-process nucleosynthesis, is also synthesized. The detec- 
tion of 7-ray lines emitted from 56 Co at unexpectedly early 
times (Matz ct al. 1988) and multi-dimensional hydrody- 
namical calculations (e.g., Fryxcll, Mullcr, & Arnctt 1991; 
Hachisu et al. 1990) suggest that the Raylcigh- Taylor in- 
stabilities occurring at each composition interface in the 
progenitor star after the passage of the blast wave have 
brought the matter in the deep region near the mass-cut 
out to the surface regions expanding at a few thousand 
km s . This mixing process finishes soon after the blast 
wave hits the stellar surface. The hydrodynamical evo- 
lution would not change the relative abundances of these 
elements inside each fluid element. Therefore it is likely 
that the spatial distributions of these three elements are 
similar throughout the explosion. With this in mind, we 
will discuss the amount of Ba synthesized in SN 1987A 
with the support of detailed observations and modelling 
for this supernova in the next section. 



0.1 



in 

>< 



0.01 - 



0.001 



0.0001 



Photspheric velocity I 
at day 20 



_i i i i_ 



_i i i i_ 



10' 



10* 



- 10 8 > 



10 7 



10 



15 



Mr Mr. 



Fig. 2. The solid curve displays the distribution of 56 Ni 
in SN 1987A inferred from X- and 7-ray observations as a 
function of enclosed mass (Kumagai et al. 1989; Shigeyama 
& Nomoto 1990). The dashed curve displays the ve- 
locity distribution of a hydrodynamical model 14E1 of 
SN 1987A (Shigeyama & Nomoto 1990). Two horizontal 
dashed lines show the expansion velocities of the photo- 
sphere inferred from observations (Menzies 1988) at day 20. 



2. SYNTHESIZED BARIUM MASS IN SN 
1987A 

The progenitor star of SN 1987A was identified as a B3 
I star Sk-69°202 (Sonncborn, Altncr, & Kirshner 1987). 
Thus the initial mass is estimated as ~20 M Q and the ra- 
dius as ^40 Rq . To estimate the stellar mass at explosion, 
one needs to consider the evolutionary path to the B3 I 
star, because this type of star undergoes mass loss due to 
the effects of stellar wind. Numerical calculations (Saio & 
Nomoto 1988) suggest that the stellar mass at explosion 
is about 16 M©. 

Multi-band photometric observations for SN 1987A 
have been used to construct so-called bolomctric light 
curves during the first 400 days after explosion (Menzies 
1988). Comparison of hydrodynamical models for this B3 
I progenitor star explosion with the bolometric light curve 
enables us to estimate the explosion energy, the amount 
of freshly synthesized 56 Ni, and the neutron star mass left 
behind (e.g., Shigeyama & Nomoto 1990). Furthermore, 
X- and 7-rays from this supernova were detected. These 
emissions originate in 56 Ni decays and the subsequent 
56 Co decays. The spectra and light curves of these emis- 
sions give us information about the distribution of 56 Ni. 
Monte-Carlo simulations of 7-ray transfer in the ejecta de- 
rived the 56 Ni distribution as a function of the enclosed 
mass (Kumagai ct al. 1989) based on a hydrodynamical 
model 14E1 (Shigeyama & Nomoto 1990) as shown in 
Figure 2. If the distribution of Ba is similar to that of 
56 Ni (i.e., the ratio Ba/ 56 Ni is uniform in the ejecta), the 
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total amount of Ba can be estimated from the mass frac- 
tion at the photosphere because both the total mass and 
distribution of 56 Ni are known. The measurement of the 
mass fraction of Ba was performed based on the spectrum 
obtained at day 20 1 and the mass fraction was found to be 
~ 2.4 x lCr 8 (Mazzali et al. 1992; Mazzali & Chugai 1995) 
. According to the above model 14E1, the photosphere at 
day 20 is located at the enclosed mass of ~12M Q (Fig. 2). 
Thus the mass of Ba is estimated to be about 6 x 10~ 6 M Q . 
A similar amount of Ba is obtained based on another hy- 
drodynamical model W10 (Pinto & Woosley 1988; Pinto 
& Woosley 1988; Pinto, Woosley, & Ensman 1988) that 
also derived the distribution of elements in the ejecta from 
comparison with observations. The, Burrows, & Bussard 
(1990) also arrived at a similar distribution to obtain the 
light curves consistent with the observed X-ray and 7-ray 
line emissions at early phases based on hydrodynamical 
model 14E1. 

3. IMPLICATIONS FROM THE DERIVED 
BARIUM YIELD 

From two different angles, we could say that the inferred 
amount of Ba ~ 6 x 10~ 6 M Q in SN 1987A is significantly 
high. First of all, this mass results in a large enhancement 
of Ba in comparison with iron (Fe) in the ejecta. The mass 
of Fe contained in SN 1987A is known with fairly good 
accuracy, because 56 Ni eventually decays to 56 Fe. As a 
result, the mass ratio of Ba to Fe in SN 1987A is found 
to be a factor of ~ 8 larger than the corresponding solar 
abundance ratio. 

Recent studies on the chemical compositions of metal- 
poor halo stars claim that these stars might inherit the 
abundance pattern of the ejecta of the preceding few 
SNe (Audouze & Silk 1995) or a single SN (Shigeyama 
& Tsujimoto 1998). To put it another way, stars are 
likely to be formed from the ISM comprising the ejecta 
of a single SN. Have we already found the very star 
that has a large enhancement of Ba as inferred in the 
ejecta of SN 1987A? The answer is yes - The stars 
are notable neutron-capture-rich giant stars CS22892-052 
and CS31082-001, whose abundance distributions of el- 
ements have been precisely determined up to Z=90-92 
(Sneden et al. 2000; Hill et al. 2001), and in which the 
neutron-capture elements are found to be overabundant: 
+0.3<[neutron-capture/Fe]<+1.8 ([X/Y] denotes the log- 
arithmic value of the number ratio of element X to element 
Y divided by the corresponding solar ratio). The Ba/Fe 
ratios in these stars are a factor of ~8-12 larger than the 
corresponding solar ratio. This coincidence of a large en- 
hancement of Ba between SN 1987A and these two stars 
might imply that these stars are descendants of SNe whose 
progenitor masses were ~20M Q . Here we compare the 
abundance distribution of some elements in CS22892-052 
with that in the ejecta of SN 1987A (Fig.l). The masses 

1 Abundance analyses of Ba were performed also for later spectra. 
However, as Mazzali & Chugai (1995) noted, the ionization by 
non-thermal electrons makes their model less plausible in deeper 
layers of the SN ejecta that become visible at later phases. 



of a few elements (Si, Ca, Ni) in SN 1987A have been 
obtained from its nebular spectra (Danzigcr et al. 1991). 
For deducing the mass of Sr, we use the same procedure 
as for Ba. The masses thus obtained have been converted 
to the abundance ratios of element X to Fe ([X/Fe]). We 
also show in this figure the abundance curve between Mg 
and Ni predicted by recent nucleosynthesis calculations for 
20M o and 13Af Q models (Umeda & Nomoto 2001). The 
abundance distribution of CS22892-052 is very similar to 
that of SN 1987A, and also to the theoretical curve of a 
20 M Q nucleosynthesis model for lighter elements. 

Secondly, the estimated mass of Ba implies that 20M Q 
SNe are the predominant sites for r-process nucleosynthe- 
sis. This can be illustrated by the following argument: 
Suppose that stars in the mass range of (20 ± AM)M & 
yield the same amount of Ba as that from SN 1987A, 
then what AM suffices to supply Ba to reach the observed 
solar r-process abundance? Since current theoretical SN 
models can predict the yields of oxygen as a function of 
the progenitor mass with fairly good accuracy (Tsujimoto 
et al. 1995; Woosley & Weaver 1995), one can estimate 
the mass ratio of these two elements ejected from all SNe 
as Ba/0~ 3 x 1O~ 7 AM/M , adopting the Salpeter ini- 
tial mass function. To explain the corresponding solar 
abundance ratio by this SN yield, we deduce AM=O.7M . 
Therefore the number of stars populating the derived mass 
range is only ~4% of all SN progenitor stars. This narrow 
mass range suggests that the abundance distributions of 
r-process elements in any star including the sun must be 
similar. This uniformity in the abundance distributions of 
r-process elements has already been pointed out by Cowan 
et al. (1999) in abundance analyses of four stars CS22892- 
052, HD115444, HD122563, and HD126238 (in the latter 
three stars, Ba is not so enhanced, i.e., [Ba/Fe]<+0.2), 
which led to the conclusion that there is only one r-process 
site in the Galaxy, at least for Z >56. Our finding is that 
the unique r-process site is (20 ± 0.7)M Q SNe, as repre- 
sented by SN 1987A. 

4. CONNECTION BETWEEN SN 1987A AND 
NEUTRON-CAPTURE-RICH STARS 

A noteworthy point is that two stars — CS22892- 
052 and CS3 1082-001 — with a large enhancement of 
Ba such as [Ba/Fe]~+1.0 have similar metallicities, i.e., 
[Fc/H] =— 3.1, —2.9, respectively. This also leads to the 
same conclusion as discussed above - that these stars in- 
herit the abundance pattern of the 20 M Q SN belonging 
to the first few generations in the Galaxy. For a given 
SN in the early galaxy, metallicities of the descendant 
stars can be estimated from the ratio of the mass of a 
certain heavy element, conventionally Fe, ejected from 
the SN, to the mass of hydrogen eventually swept up by 
the explosion (Shigeyama & Tsujimoto 1998). Neutron- 
capture-rich stars must be the first few generation stars, 
because in later generation stars, such a large enhance- 
ment of r-proccss elements would be reduced due to a 
gradual increase in the contribution from the ISM with 
low [r-process/Fe] ratios to the stellar abundance. Since 
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each SN with a different progenitor mass yields different 
amounts of heavy elements, we obtain the stellar metallic- 
ity (or the Fe/H ratio) as a function of the progenitor mass 
of a SN from which the star was born. Taking a ratio of 
the Fe yield ~0.07M Q to the swept-up mass of hydrogen 
~ 6.5 x 10 4 Mq for a 20 M SN, we obtain [Fe/H]- -3. 
The metallicities of CS22892-052 and CS31082-001 indeed 
suggest that their progenitor masses were ~20 Mq. 

Once the progenitor mass for these stars has been de- 
termined, we can estimate the mass of Ba ejected by a 20 
Mq SN, combining the observed [Ba/Fc] ratios in these 
stars with the Fe yield from the SN. In fact, it has been 
already found (Tsujimoto et al. 2000) that CS22892-052 
was born from a 20 M Q SN that had produced 8 x 10~ 6 
Mq of Ba. Thus there is a fairly good agreement in the Ba 
mass ejected from a 20 Af© SN independently predicted 
by both SN 1987A and CS22892-052. It is a natural result 
of their similar elemental abundance distributions. 

Finally we need to examine Sr. The overabundance of 
Sr in comparison with Fe in SN 1987A is relatively small, 
compared with the prediction by theoretical r-process line 
(Kappler, Beer, & Wisshak 1989) in Figure 1. This might 
imply that Sr has other production sites besides 20M Q 
SNe. This consideration is compatible with the observed 
breakdown of the concordance between the solar abun- 
dances and the CS22892-052 abundances for the lighter 
r-process elements (Sneden et al. 2000). One possible 
candidate for other sites may be a wider mass range for 
the production site of Sr, e.g., (20±3)M Q . This range 
is obtained from the same procedure for estimating the 
value of AM as was performed for Ba. The obtained AM 
for Sr approximately corresponds to the metallicity range 
—3.3 [Fe/H] —2.7. It implies that among stars in this 
metallicity range, there exist stars in which Sr is much en- 
hanced compared with Ba. As a matter of fact, HD122563 
([Fe/H] =-2. 74) has a very large enhancement of Sr such 
as [Sr/Ba]=+l.l (Wcstin et al. 2000). 

5. CONCLUSIONS 

We have shown that r-process nucleosynthesis is cer- 
tainly evident in SN 1987A. Our study clarifies that the 
observed overabundance of Ba in the outer layer of the 
ejecta of SN 1987A is part of the r-process production 
synthesized during explosion in the deepest layers of the 
ejecta, though this overabundance was ascribed to the re- 
sult of the s-process operating in the core He burning 
stage. We presented three strong arguments for this. 

First, we should pay attention to the observed fact that 
the overabundance of Sr was also found at the same time 
in SN 1987A and that the abundance of Sr is less enhanced 
than that of Ba, which makes it unlikely that these ele- 
ments are of s-process origin. In addition, the derived 
Ba/Sr ratio in SN 1987A is rather similar to those ob- 
served from metal-poor stars in which Ba has been shown 
to be of r-process origin from the observed Ba/Eu ratios. 

Secondly, the absence of an overabundance of Ba in 
other SNe II-P contradicts the prediction from theoretical 
models that the s-process operates irrespective of He-core 



masses. Alternatively, this fact implies that the r-proccss 
occurs in SNe whose progenitor masses are limited to a 
narrow range. Such a unique production site for r-proccss 
is fully consistent with the uniformity in the abundance 
distributions of r-process elements among several stars in- 
cluding the sun. 

Finally, the total amount of Ba inferred from the Ba 
abundance in the outer layer of the ejecta from the stand- 
point of the r-process origin reveals that 20 Mq SNe, one 
of which is SN 1987A, are the prominent production sites 
for r-process elements. This is consistent with the sec- 
ond point above, i.e., no overabundances of Ba in type 
II-P SNe other than SN 1987A and the uniformity in the 
stellar abundance distributions of r-process elements. 

We have further shown that the elemental abun- 
dance distributions of neutron-capture-rich giant stars — 
CS22892-052 and CS31082-001 - are very similar to that 
of SN 1987A, which suggests that these stars are born 
from the ejecta of 20 M Q SNe like SN 1987A in the early 
Galaxy. The high abundances of Ba observed in these 
stars also imply that the major production sites for r- 
process are —20 Mq SNe. 

This work has been partly supported by COE research 
(07CE2002) and a Grant-in-Aid for Scientific Research 
(12640242) of the Ministry of Education, Science, Culture, 
and Sports in Japan. 
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